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Abstract 

For an m-dimensional multivariate extreme value distribution there exist 2™ — f exponent 
measures which are linked and completely characterise the dependence of the distribution 
and all of its lower dimensional margins. In this paper we generalise the inequalities of 



Schlather and Tawn ( 2002 ) for the sets of extremal coefficients and construct bounds that 
higher order exponent measures need to satisfy to be consistent with lower order exponent 
measures. Subsequently we construct nonparametric estimators of the exponent measures 
which impose, through a likelihood-based procedure, the new dependence constraints and 
provide an improvement on the unconstrained estimators. 

Keywords: max-stable distributions; multivariate extremes; exponent measure; inequalities; 
constrained estimators 

1 Introduction 

Max-stable distributions arise naturally from the study of limiting distributions of appropriately 
scaled componentwise maxima of independent and identically distributed random variables. 
Here and throughout the vector algebra is to be interpreted as componentwise. A vector random 
variable Y = with unit Frechet margins, i.e., Gi{y) := ¥{Yi < y) = exp(— 1/y), 

y > 0, i & Mm = {1, . • • ,m}, is called max-stable if its distribution function is max-stable, i.e., 
if 

GMmiVMm) ■■=^{Y < yMra) = Bxp <^ - / max ( — ) dH{wi, Wm) \ , (1) 



where yM„ = ivi, Vm) e Sm = {{wi, . . . , Wm) G K^J" : YT=i = l} is the (m - 1)- 
dimensional unit simplex and H is an arbitrary finite measure that satisfies 



/ WidH{wi, . . . , Wm) = 1 for any i € M„ 

'5'm 



The last condition is necessary for to have unit Frechet margins and representation ([T]) is 



due to Pickands (1981). There is no loss of generality in assuming unit Frechet margins since 



our focus is placed on the dependence structure of max-stable distributions, i.e., we are inter- 



ested in the copula function (Nelsen, 1999) which is invariant to strictly monotone marginal 



transformations and in practice we can standardise random variables to unit Frechet margins. 



The dependence properties of max-stable distributions have received attention in the multivari- 



ate extreme value literature. Dating back to Sibuya (1960) and Tiago de Oliveira (1962/63), 
it has been known that max-stable distributions are necessarily positively quadrant dependent, 



I.e., 



(2) 



1=1 



which implies that no pair of random variables can be negatively dependent. Additionally, 



max-stable distributions satisfy even stronger forms of dependence. Marshall and Olkin ( 1983 ) 



show that Cov{g{Y), h(Y)} > for every pair of non-decreasing real functions g and h on M'", 
i.e., they are associated. For a review of the dependence properties of max-stable distributions 



we refer the reader to Beirlant et al. (2004) and the references therein. 



Although all of the aforementioned properties exhibit characteristics for the dependence struc- 
ture of the class of max-stable distributions, they are far too general to be either tested or 
implemented in practice. In this paper, we introduce additional constraints for the dependence 
structure that can be incorporated, through a likelihood-based procedure, into the estimation 
of max-stable distributions from observed componentwise maxima. The new constraints are in 



essence the generalisation of the Schlather and Tawn ( 2002 2003 ) inequalities for the extremal 



coefficients which correspond to the dependence properties of max-stable distributions for the 
special case of GMmiu^ • • • i y)) y > 0. As such, our notation and strategy are infiuenced by the 
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work of Schlather and Tawn (2002, 2003[ ). The new inequalities presented in this paper are 
related to the general case of GMm{yMm)i UMm ^ '^+- 



In Section [2] we introduce the class of max-stable distributions along with the Schlather and 



Tawn ( 2002[ ) inequalities for the extremal coefficients. Subsequently, we present the general 
result of the paper that gives rise to inequalities for the exponent measures. In Section |3] 



we consider the Hall and Tajvidi (2000) nonparametric estimator for the exponent measure 



and extend it, through a likelihood-based procedure, to satisfy the new inequalities. Finally, in 
Section [4] a simulation study is conducted to assess the performance of the constrained estimator. 



2 Dependence Properties 
2.1 Background 

The class of max-stable distributions arises naturally from the study of appropriately scaled 
component-wise maxima of random variables. Consider a set of independent and identically 
distributed random vectors = {Xl, . . . , Xm), j = 1, . . . , n, with unit Frechet margins. Under 



weak conditions (Resnick, 1987) it follows that 

GMr^iVMj, VMrr^eR^. (3) 



The distribution function Ga/^ can be completely characterised by the following representations 

Vm„ {yMm ) = - log GMrn (y) = / max ( — ) dH{wi, . . . , Wm) , (4) 

Jsm V Vi J 



(5) 



where the function Vm,„ is known as the exponent measure of the multivariate extreme value 
distribution G Mm and Am^ ; called the Pickands ' dependence function, is a convex function that 
satisfies ma,x{wi, . . . , Wm} < ^Mmiwi, ■ ■ ■ , Wm) < 1, {wi, . . . , Wm) £ Sm- This condition implies 
that ^Mm(ej) = 1, J G Mm-, where Cj is the j-th unit vector in W^. 
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Let Cm = 2*^" \ {0} and denote also by i/b = {m '■ i S B} for B G Cm- Then we can define 
2"^ — 1 exponent measures for an m-dimensional max-stable random vector Y , where each one 
characterises completely the distribution function of a marginal random variable Yb of Y , i.e., 



VsiyB) = -\og{¥{YB < ys)] = - 




The set of exponent measures {Vb '■ B ^ Cm} describes completely the dependence struc- 
ture of a max-stable distribution given by equation ([T]) and all of its lower dimensional mar- 
gins. It is also trivial to see that with each exponent measure Vb there is an associated 
Pickands' dependence function Ab- Additionally, Vb, B £ Cm, is homogeneous of order — 1, 
i.e., VB{y,...,y) = y-^VB{l,...,l),y>0. 



The importance of the homogeneity property is mostly illustrated through one widely used 
measure of extremal dependence for the variables indexed by a set B £ Cm- More specifically, 
the quantity defined by 

eB = VB{l,...,l)= [ msiXWidH{wi,...,Wm), 1<0b<\B\, (6) 

describes the effective number of independent variables in the set B and arises naturally from 
the distribution of the maximum of all the variables indexed by the set B, i.e., 

Pjmaxy, <y| =P{y, <yfs^ ^ > 0. (7) 

The measure is termed the extremal coefficient and complete dependence and independence 
corresponds to 9b = 1 and 6b = \B\ respectively. Also, from expression ([T]) it follows trivially 
that 9b = 1 for any B G Cm with \B\ = 1. Due to its simple interpretation, the set of extremal 



coefficients {9b : B £ C} has been used as a dependence measure in various applications (Tawn 



1990 Schlather and Tawn 2003) 



2.2 Schlather and Tawn (2002, 2003) inequalities for the extremal coefficients 

Schlather and Tawn (2002, 2003) constructed bounds for the set of extremal coefficients {9b '■ 



B € Cm} of max-stable distributions that characterise the dependence structure for the special 



case ofGMmiy^ ■ ■ ■ ,y), U > 0- Here we use the terminology of Schlather and Tawn (2002) and for 
non-empty distinct subsets Bi, . . . , Bg of Mm, s € N, we refer to the set of extremal coefficients 
{9bi, ■ ■ ■ ,(^Bs} as complete and consistent if s = 2"^ — 1 and 0^. is given by expression ([6]), 
respectively. Their main result is given in the following theorem. 



Theorem 1 (Schlather and Tawn (2002) Corollary 5). A complete set of extremal coefficients 



{9b '■ B ^ Cm}, where Mm is a finite set of indices, is consistent if and only if 

J2 > 0, for all LC^ Cm. (8) 

BeC,n,BDMm\L 

Theorem [T] yields bounds that higher order extremal coefficients need to satisfy to be consistent 
with lower order extremal coefficients. For example consider the inequalities ([s]) for the cases 
m = 2 and m = 3 and let for ease of notation O^^ jy and d{ij^k} be 6ij and 9ijk, for i,j, k G Mm 
and i ^ j ^ k. These are respectively 

1 < 9i2, 9i3, 023 < 2 
and 

max {012, 013, 023, O12 + ^13 + O23 - 3} < 6*123 < min {0i2 + 0i3 - 1, O12 + ^23 - 1, 6*13 + O23 - 1} • 

The first set of inequalities represents the well known bounds of the extremal coefficients that 
come from the positively quadrant dependence property ^ of max-stable distributions. How- 
ever, the second set of inequalities gives tighter bounds for the higher order extremal coefficient 
^123. This can be seen easily since the combined inequalities for the cases m = 2 and m = 3 
reduce to 1 < ^123 < 3. 

2.3 Inequalities for the exponent measures of max-stable distributions 

It transpires that similar inequalities as with those in expression ([s]) can be obtained for the ex- 
ponent measures {Vb ■ B £ Cm} of max-stable distributions. Analogously with the terminology 



for the extremal coefficients in Section 2.2 we introduce the following definition 



Definition 1. Let s be an integer, for i = 1, . . . , s Bi are distinct non-empty subsets of M„ 
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{1, . . . ,m} and yu^ = {vi, ■ ■ ■ ,ym) S M+. An ensemble {Vs^iyBi), Vs^CysJ} of exponent 
measures, where = {yj ■ j G -Bj}, is called consistent if 



Vb,(2/bJ = / maxl^] dH{wi,...,Wm), 

J Sr. je^' V yj J 



for i = 1, . . . ,s and H is an arbitrary finite measure that satisfies Jg WidH{wi, . . . , Wm) = 1 
for any i £ Mm- 

If s = 2™ — 1 then the set of exponent measures is caUed complete. The foUowing theorem pro- 
vides a new representation of the exponent measures of multivariate extreme- value distributions 
in terms of non-negative and uniquely defined real functions. 

Theorem 2. Let {Vb ■ B € Cm} be a complete and consistent set of exponent measures. Then, 
there exist 2™ — 1 non-negative functions dL ■ — )• M+, L G Cm, such that, for any B G Cm 

VB{yB)= dLiyMj, (9) 

and the functions di are uniquely given by 

dL {yuj = Yl (-l)!'"'"''!"^' / max ('^) dH{w,, . . .,Wm). (10) 
^-^ Q jeB \ Vj I 

Proof 



The proof of equation Q of Theorem [2] follows along the lines of Schlather and Tawn ( 2002 ) 



proof of Theorem 5 for the simpler case of the extremal coefficients by replacing the constants 



ot\{n) of Deheuvels (1983) representation of max-stable distributions with a\{n)/yi, i € Mm, 
k £ Ij. Equation (10) is the Mobius inversion of equation ([9]).n 

The characterisation of a consistent set of exponent measures is obtained from the following 
corollary. 

Corollary 1. A complete set of exponent measures {Vb '■ B G Cm} is consistent if and only if 
V (-l)l^nL|+i f fwA ^ . . . , > 0, (11) 

B&Cm,BDMm\L -J^m-I \ / 

for all yUm ^ ^+ (^^d L G Cm- 



3 Inference 



3.1 The Hall and Tajvidi (2000) estimator of the exponent measure 



The fundamental premise in all statistical extreme value modelling is that the observed extremes 
of a stochastic process are well modelled by the limiting theoretical extreme- value distributions. 
Let for example = {X^^ . . . j = I, . . . ,N, be a set of independent and identically 

distributed m-dimensional random vectors with unit Frechet margins. Here and throughout we 
assume that the normalised componentwise block maxima 



V 

r=ij-l)d+l 



X' 



-, j = l,...,n, 



where nd = N, follow exactly the law Gm,^ of the limiting expression Q. 



Let now wb G 5*1^1 ieB'^B,i - l|, -B G Cm, and define Z^^ - wbY^, for 

,n. It then follows that the cumulative distribution function of maxjg^ is Frechet 



J 



1, 



with scale parameter equal to the Pickands' dependence function Ab{wb) of Gb, i-e., 



maxZ] < y 



exp < ^ } , y>0. 



A natural consistent estimator of Ab then is the Hall and Tajvidi (2000) corrected version of 



Pickands' estimator (Pickands, 1981) which maximises the likelihood 



{Ab{wb)} = nlog{AB{wB)} - 2Y,^ogWi - Ab{wb)^—j, 

i=i i=i 



(12) 



Hall and Tajvidi 



(2000) 



where = maxjgB ^WB,iY/ Y!j=i{l / ) / n |, j = 1, . . . , n, is the 

correction which ensures that maxw^ < -^b{wb)-, for all wb G S\b\i as well as Asiej) = 1, 
for any j G M^, where Cj is the j-th unit vector in M™. The maximum likelihood estimator is 
given by Ab{wb) = l^"-*^ Sj=i(l/^B)| which is subsequently corrected by 

^b{wb) = min ^Ab{wb), i| 



to satisfy Ab{wb) < 1, for all wb G S\b\- On combining the estimator Ab with equation (^5^, 



the following consistent estimator of the exponent measure Vb is obtained, 



VB G 



,\B\ 



B GCr, 



(13) 



Other types of estimators exist in the literature such as the non-parametric estimators proposed 



by Deheuvels (1991) and Caperaa et al. (1997) for the bivariate case. Zhang et al. (2008) gives 



a detailed overview of the existing estimators and extends them to the multivariate case. In 



this paper though we use the Hall and Tajvidi ( 2000 ) estimator since it arises as the maximum 



of a log-likelihood function based on which the new inequalities of Section 2.3 can be imposed. 



3.2 Constrained estimators 

It transpires that the aforementioned nonparametric estimators of the exponent measures do 
not necessarily ensure that the resulting estimated set of exponent measures satisfy inequali- 



ties ( 11 ). The focus here is placed on incorporating these additional constraints in the estimation 



procedure so that the resulting complete set of estimated exponent measures VsivB), B G Cm, 
is consistent in the sense of Definition [T] for fixed G M™. To incorporate the inequali- 



ties we construct similarly with Schlather and Tawn (2003) a joint log-likelihood function £ of 



{Ab]B G Cm} by falsely assuming independence between the observations for all different B to 
give the pseudo-log-likelihood 



AB{^^^):B^Cm])= E 



(14) 



The maximum pseudo-likelihood estimators are consistent (Liang and Self, 1996) and the con- 



strained estimators are obtained by maximising the pseudo-log-likelihood ( 14 ) subject to 



Bee, 



E (-i)'^""'^^(evi/.Ub 

m,BDMm\L {i&B J 



B 



> 0, for all L G C„ 



and 



At 



< 1, for all 5 G C„ 



The resulting constrained estimators are denoted by A'^ which in turn yield the estimators 



as in equation (13). The joint estimation of the exponent measures ensures that all estima- 
tors are self-consistent. Note that the resulting estimates of lower order exponent measures 
are affected by higher order measures, i.e., estimates of Vbo(?/Bo) affected by estimates of 
Vb^ (y_Bi)i where Bq C Bi. The major benefit of this feature is that this guarantees the existence 
of higher order measures which are self-consistent with the lower order measures. 

An alternative way of obtaining a set of estimated exponent measures is via sequential es- 
timation, i.e., the lower order exponent measures are estimated firstly and then are used as 
constraints in the estimation of the higher order exponent measures, see also |Schlather and| 



Tawn (2003). Although this method is faster than the joint optimization problem described by 



equation (14), it does not have the desirable feature described above. 



4 Simulation Study 
4.1 Design 



We illustrate the impact of constraining the Hall and Tajvidi (2000) estimators to satisfy the 



new inequalities (11) over the unconstrained estimators of the set of exponent measures using 
simulated data from a 3-dimensional max-stable distribution, i.e., the extreme value logistic 
distribution with dependence parameter a G (0, 1] and set of exponent measures given by 

|v^B(yB)= (^l^y,"'^"^ ■■ BeC^y vb^R^^^. (15) 

The values a = 1 and a = 0, taken as a — )• 0, correspond to independence and complete depen- 
dence, respectively. 

All comparisons are based on the root mean square error (RMSE) performance of the exponent 
measure estimators for a range of dependence parameters a and a cube grid of values, say 
L'^ C M.^^ for t/Mg- Specifically, the values chosen for the dependence parameter and the sample 
size are a S {0.2, 0.5, 0.8} and n = 50, respectively. Results from larger sample sizes are not 
reported in the paper since they are unrealistic for applications and also, the efficiency of the 



estimators Vb and Vn is similar, a fact that comes from the consistency property of the Hall 
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and Tajvidi (2000) estimator. The set L was chosen to be the discrete set {xp^, . . . , Xp^} with Xp 
denoting the p-th quantile of the unit Frechet distribution. We chose pi = 0.05, p7 = 0.95 and 
step size pj —pj-i = 0.15. The Monte Carlo size used to compute estimates of the RMSE is 500. 

To obtain an aggregated measure of performance, we also report the Monte Carlo estimates of 
the integrated square deviation of the estimators from the theoretical function, i.e., 

fB= [ \vB{yB)-VB{yB)Y dvB, for ah B e C3, (16) 

Jc{h\B\) ^ J 

where C(LI^I) is the smallest |i?|-hypercube that contains the set LI^L The integral in expres- 



sion (16) is approximated in each Monte Carlo iteration by the quadrature mid-point numerical 
integration technique on the grid L'^' and the measure is defined analogously by replacing 
Vb in expression (16) with V^. 



4.2 Results 

Figure [1] shows the histograms of the ratio of RMSEs between and Vb, B £ C3, for all grid 
points in M.^ for the extreme value logistic distribution with a = 0.2,0.5 and 0.8. The figures 
indicate either similar or better performance of the constrained estimators. 
In particular, for the a = 0.8 case, the constrained estimators are more efficient than the un- 
constrained estimators especially for the higher order exponent measure V123 and improvement 
in RMSE, although lower in magnitude, can be also seen in the bivariate exponent measures 
Vb, B G C3\M3. Also, the percentage of Monte Carlo samples where the constrained estimates 



changed with respect to the Hall and Tajvidi (2000) estimates is 62%. Regarding the a = 0.5 



case, we found better performance of the constrained estimators for V123, although lower in 
magnitude than the a = 0.8 case, and similar performance for the bivariate exponent measures. 
This feature is also supported by the smaller percentage of change in estimates which is 30%. 
For the case of strong dependence, i.e., a = 0.2, the percentage of change in estimates is very low 
and equal to 6% which results in similar efficiency of the estimators for all exponent measures 
as is also shown from Figure [T| 

Table [T] shows the Monte Carlo estimates of the integrated square deviation of the estimators 
from the theoretical function. For the case of strong dependence there is no practical benefit 
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Figure 1: Histograms of the ratio of Monte Carlo estimates of RMSEs between the constrained 
and unconstrained estimators of the exponent measures F12, V13, V23 and F123 (top to bottom) 
for the extreme value logistic distribution with a = 0.2 (left) a = 0.5 (centre) and a = 0.8 
(right). 



Table 1: Monte Carlo estimates of Tg and T^, B & C3, for the extreme value logistic case with 
a = 0.2,0.5 and 0.8. 





a = 


0.2 


a = 


0.5 


a = 


0.8 


B 


Tb 




Tb 


T'b 


Tb 


T'b 


{1,2} 


0.02 


0.02 


0.36 


0.35 


0.76 


0.72 


{1,3} 


0.02 


0.02 


0.34 


0.34 


0.76 


0.74 


{2,3} 


0.02 


0.02 


0.37 


0.36 


0.74 


0.72 


{1,2,3} 


0.66 


0.66 


11.15 


10.63 


26.80 


22.14 
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of over Vb- However, in all other cases the constrained estimators are more efficient than 
the unconstrained estimators. This shows that not only does the imposition of the constraints 
improve the performance of the estimators for the higher order exponent measures, but so does 
for the bivariate level of dependence. 

To conclude, the performance of the estimators Vb and is similar as the dependence increases 
and becomes identical in the limiting case of a — )• 0. This feature is explained by the increase 



in performance of the Hall and Tajvidi (2000) estimators Vb as dependence increases which 



yields a consistent set of estimated exponent measures. Overall, we found the imposition of 
the new constraints to be beneficial for the simplest max-stable distribution, i.e., the extreme 
value logistic, and superior in efficiency, especially for the case of moderate or weak dependence. 
The largest improvement is observed for higher order exponent measures which is promising for 
implementations in higher than 3 dimensions. 
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